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Abstract 



We study the Laughlin wave function on the cy Under. We find it only 
describes an incompressible fluid when the two lengths of the cylinder are 
comparable. As the radius is made smaller at fixed area, we observe a con- 
tinuous transition to the charge density wave Tao-Thouless state. We also 
present some exact properties of the wave function in its polynomial form. 
We then study the edge excitations of the quantum Hall incompressible fluid 
modeled by the Laughlin wave function. The exponent describing the fluctu- 
ation of the edge predicted by recent theories is shown to be identical with 
numerical calculations. In particular, for = 1/3, we obtain the occupation 
amplitudes of edge state n{k) for 4-10 electron size systems. When plotted as 
a function of the scaled wave vector they become essentially free of finite-size 
effects. The resulting curve obtains a very good agreement with the appropri- 
ate infinite-size Calogero-Sutherland model occupation numbers. Finally, we 
numerically obtain n{k) of the edge excitations for some pairing states which 
may be relevant to the u = 5/2 incompressible Hall state. 
PACS numbers 73.40.Kp, 73.20.Dx, 73.50.Jt 
June 1994, v2.6 
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It is now well-established that the quantum Hall effect (QHE) in the spin-polarized 
two-dimensional electron liquid at Landau-level filling-factor u = 1/q results from an incom- 
pressible correlated electron state which is very well described by the simple model wave 
functions introduced by Laughlin This has been verified by extensive numerical studies 
of systems of a finite number of electrons which have been carried out in two popular 
geometries: spherical [Q] and periodic, or toroidal 

These geometries are convenient for obtaining bulk properties by extrapolation to the 
thermodynamic limit as they do not introduce edges. In the most symmetric spherical 
geometry, the geometry is fully specified by the requirements of translational and rotational 
symmetry; in the periodic case, rotational symmetry is absent, and while the area of the 
elementary cell is fixed by the integrally-quantized number of magnetic flux quanta that pass 
through the surface, the Bravais lattice of translations in which physical (gauge- invariant) 
quantities are formulated can be continuously varied. When the shortest of these translations 
has a length comparable to or shorter than a magnetic length (or than the mean bulk 
interparticle spacing) one may expect that the results of the finite-size calculation are no 
longer representative of the bulk physics. 

There are other geometries in which the QHE has been discussed; in particular, cylin- 
drical geometry appears well suited for studies of edge states associated with the boundaries 
of a region of incompressible fluid. In this geometry, there are two characteristic lengths in 
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the A^-electron system: one is the periodic repeat distance (circumference of the "cyhnder" ) 
which plays the role of the shortest translation in the periodic system; the other is the width 
of the finite strip of incompressible fluid that results when the electrons are confined in the 
well of a potential that varies in the direction parallel to the axis of the cylinder. When 
this length becomes small, there are interactions between the two edges of the fluid, and 
calculations in this limit (unlike the case of small periodic repeat distance) are presumably 
relevant to the physics of the QHE in narrow channels. 

In this paper we describe the properties of the Laughlin wave function on the cylinder 
as the circumference of the cylinder is varied at fixed particle number (i.e., fixed surface 
area of the incompressible fluid) and study its' edge excitations. As usual, we only consider 
electrons restricted to states within the lowest Landau level. 

The Laughlin wave function on the cylinder 0| for z/ = 1/m is given by: 

2 

-^L = \{{e— -e—) n^^, (1) 

i<j i 

where 7 = ^/-R, £ = (h/eB)^ is the magnetic length, and R is the radius of the cylinder. The 
complex coordinate z is x + iy, where x is the periodic coordinate along the circumference 
of the cylinder and y is along the axis. 

The periodicity in the x-direction means that the the allowed pseudo-momentum kx of 
each particle is quantized: = n/R, where the integer n labels the conventional basis 
of single-particle "Landau gauge" states. We allow the range of y to become infinite, but 
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"compress" the electrons by restricting them to the Landau-gauge states with n in the range 
0, 1, 2, . . . , N^. Just as in the spherical geometry, the Landau level degeneracy is equal to 
Nl = + 1, which is proportional to the area available for the guiding center motion. 

It is important to note that the study of many-particle states on the cylinder is different 
from other geometries because of a similar degeneracy if the ground state is an incompressible 
fluid. In this case, the fluid will determine its own density and need not occupy the entire 
surface available to it, but will form a "ribbon" on the cylinder. The displacement of this 
ribbon as a whole will cause an overall degeneracy. For finite systems interactions with the 
edges will lift these into quasi-degeneracies. As a result, care needs to be exercised in sorting 
out the internal excitations from the translational modes. 

We construct the Laughlin wave function numerically as the only zero energy state of 
the hard-core potential with pseudo-potential parameters 7^ for m = 0, 1, . . . , z/~^ — 1, 
where Vm are the pair energies in the state of relative angular momentum m. 

The Hamiltonian can be obtained by a generalization of the pseudo-potential formulation 
to the cylinder as follows: 

H=^22z2 dqVmLm{q'r + rn')e 2 e 2 e ^ e*^^^" ^^^e ^ . (2) 

The parameter 7 serves as an "aspect ratio" (we fix the area to preserve the Landau level 
degeneracy), and it proves to be the crucial parameter in determining whether the ground 
state is an incompressible fluid or not. For finite systems, the physical properties of 
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depend strongly on 7. For 7 > 0.7, becomes different from an incompressible uniform 
fluid and begins to approach a charge density wave-like state, which for larger 7 (~ 1.2) is 
essentially described by the Tao-Thouless 0] (TT) state. 

To investigate the 7 dependence, we note that in occupation space has the general 
form: 

^L=T, A{n,, n2, . . . , un) Y[ {e'^'T 11 (3) 

{nj i i 

-i: /'"'i:;));""' n(e^)n^..w w 

. A(ni,n2, . . . , n^r) are amplitudes independent of 7, 0(7) is an overall normalization factor 
independent of the occupation numbers, and "^niz) are the appropriate single particle wave 
functions. An algorithm for calculating these amplitudes for small system sizes has been 
given previously [Q, but this does not appear to provide any substantial advantage to the 
alternative of diagonalizing H directly. Accordingly, we calculate numerically from H 
for 7 = and then construct it for 7 7^ from Eq. 4. The charge densities shown in Fig. 
1 and 2 were obtained in this manner. The parameter a in the figures is the aspect ratio: 
a = L/[2t[R) = (7^ X Nl)/2tt, where L is the length of the cylinder. The square points 
in Figs. 1 and 2 are the average occupation of each orbital positioned at the center of the 
Gaussian envelope of the orbital, which are separated by Ay = 7. These figures capture the 
transition to the charge density wave state as 7 is increased. 

The extremely squeezed cylinder (hoop geometry) is realized at 7 = 0; this is an inter- 
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esting limit where the spatial dimensions have been reduced to essentially one in that the 
Gaussian envelopes are always centered at y = 0. The density however does have a finite 
width in the squeezed direction (see Fig. 3). We have found several interesting properties 
of the Laughlin's wave function for 7 = which we now address. In this limit, assumes 
the familiar form: 

i<j i 

Clearly in the expansion of this wave function all occupation amplitudes are integers. In 
particular, expanding the Jastrow factor we obtain: 

N 

]li^^-^,r= E C^(mi,m2,...,m^)n^r- (6) 

i<j mi,m2,---,mn i=l 

In the appendix, we show that: 

• (a) C{{{mpi}) — ((— l)^)"*C(({mj}), P is a permutation on N objects. 

• (b) C({m°}) = 1, where {m°} = {m{i - 1)}. 

• (c) C({mj}) = mx integer, (integer can be zero) if {rrii} can be obtained from a 
permutation of {m°} by a succession of "squeezing" operations (defined below). 

• (d) C({mi}) = if {rui} is not a permutation of {m°} or a configuration obtained by 
successive "squeezing" of {m°}. 

The squeezing operation (without reordering) {rrii} is defined for some pair 

j ^ k with nij > nik and a positive integer n by: 



m ■ = rui if, i ^ j,k 

> where < n < rrij — m^. 



m'j = rrij — n 



Note that: 



^(mj - rrij)^ > '^{m[ - m'^f = ^{rrii - rrijY + 2nN{n - {rrij - rrik)). 

i<j i<j i<j 

Since every occupation state with non-zero amphtude can be squeezed from the TT state 
and, since within constants, J2i is related to I]i<j(^i ~%)^; it can be seen from Eq. 4 that 
the above inequahty imphes that the TT state will dominate the wave function as 7 — > oo. 
In practice this condition is realized for 7 > 1, see Fig. 2. 

We have also empirically verified the following for up to A?" = 8 electron systems, the 
proof of which we leave as an open problem: 

• (1) The filled level droplet has the highest amplitude. 

• (2) ^(filled level droplet) = ^''^^^y^^' 



(3) The total normalization is ||^| 



The first four items (a-d) are also vahd for bosons except of course the wave function is 
now symmetric. The last three items (1-3) are however only valid for fermions. 

It should be noted here that in the squeezed limit the appropriate 1-d Hamiltonian for 
which the u — 1/3 Laughlin state is an exact zero-energy ground state is: 

H{xi,X2,---,XN)^-Y.-Q^S{Xi-Xj) (7) 

i<j i 



This result also noticed by Wen is an obvious limit of the full two-dimensional hard core 
potential [gl. 

Before we study edge states we would like to make a few comments. The period of the 
oscillation of the density near the edges in the fluid phase is determined by the bulk density- 
density response function which is dominated by a single mode namely the magneto-roton 
||1CI|| . As such, it is similar to the impurity screening (or the lack of it) by the incompressible 
fluid. The period at the edge is comparable to that of the density oscillations near an 



impurity [111]. In the limit of large aspect ratios, the system makes a transition to a "stripped 
Wigner Crystal" phase where the number of peaks in the density is equal to number of 
electrons. This same limit, but for the Coulomb interactions, has been studied by Chui 



12 in connection with narrow channels. However, he used toroidal boundary conditions. 



It is not surprising to see gapless modes |T3| (as we see) in this limit. In the opposite 



limit, e.g. small 7, the gap survives at least for the hard-core repulsion and may be the 
more appropriate limit for narrow channel geometries. Whether or not QHE can be ruled 
out for very narrow channels has not yet been addressed in our study and needs further 
investigation. 

We next turn to the discussion of edge excitations of the Laughlin droplet (or ribbon 
in our case). Again, the most interesting limit turns out to be the completely squeezed 
cylinder (7 = 0). In this limit there is no bulk in the usual sense and only edge excitations 



remain. Wen has used current algebra techniques to study edge states with the same 



results as here. One of us [T^ has interpreted the edge states as "generalized Fermi surface 
singularities" , obeying a generalized Luttinger [|l^] theorem. Here we follow that formalism. 

Consider first a free Fermi gas with a momentum space energy surface E{k) having the 
topology shown in Fig. 4. In terms of the Fermi points A;^, we have the Luttinger theorem 
for the total number of particles and the total momentum of the system P: 



iV = A^A..A:} (8) 

i 

p=^j:^Mk}r (9) 



where Au = n{kf — e) — n{kf + e) is the jump across the Fermi point, and L is the length 
of the system. One can generalize these expressions to contain local deviations of the Fermi 
points 6kj{x) from their uniform value ky. 



^p{x) =—Y.^^^^{x) (10) 

An(x) = ^ E ^l^^k}Sk}{x) + ^Mm^W (11) 



where, 



k){x) = k) + ^ (12) 
p{x) = J (13) 
U{x) = J. (14) 
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The local description is valid only if for any pairs of indices i,j the condition (fcj — fcj.)^ ^ 1 
is satisfied. Here ^ is an inverse momentum cut-off measured from the Fermi points. 

Similarly the Hamiltonian for edge deformation can be expressed in terms of the local 
fluctuations of the Fermi wave vector: 

H^'^ = ^J dxY.v)^v,{5k){x)f (15) 

i 

where f j's are the Fermi velocities. It is possible to include interactions terms between these 
fluctuations as in Fermi Liquid theory but such terms are not relevant to this work and will 
not be pursued here. We next proceed to quantization. 

It is well known that for 1-d Fermi gas |T^, the Fourier components of the density 



operators (Tomonaga bosons) 



P\^) = \ E e^'V; (16) 

satisfy the Kac-Moody algebra: 

K Pi'l = ^^Aj^S.+i'fi, < 1. (17) 

Furthermore, Bosonizations of Fermi fields in 1-d give: 

^l(x) = A,e* (18) 
= Aic"^^^''^ (19) 

where Aj, to be given shortly, are operators which make \&'s anticommute, and <^j(x) is the 

canonical Bose field. From these and Eq. 12 we obtain: 
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-4^^ (21, 



which is the "chiral constraint". Next, from the above Eqs., we obtain the commutation 
rules (CR): 

[p,ix),p,{x')] = S,,^^{S{x-x')) (22) 



[P^(^)'^]=^M^(^(^-^')) (23) 



{x),p,ix')]=tS,Jix-x') (24) 



(x), = 2mAu,6,,,6{x - x') (25) 



1J[ 

[(l)i{x),(l)j{x')] = —6ijSgn{x - x') (26) 



Thus 0, and p (or the derivative of 0) are conjugate fields. Using These CR's and 

[N,¥,{x)] = ¥,{x), (27) 

it can be seen that: 

<i!l{x)^i{x') = ^i(a;')^I(a;)e-['^'(^)''^'(^')l. (28) 

For the fields to obey Fermi statistics we must require that ^ be an odd integer ±m. We 
also note that the choice 
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A, = e^^^^^"('-^'^'^^- (29) 

makes the Fermi fields for different sectors {i 7^ j) anticommute as well. 

The usual results for the 1-d Fermi gas are recovered when Au = 1. On the other hand, 
applying this formalism to the quantum Hall effect we obtain Au = 1/m. The CR's for the 
Bose fields in these two cases are therefore related by: 

0(a;')]QHE = ~['^(^)' <^(^')]Fermi (30) 



Following standard treatment of the Luttinger model ^^7j , except here one needs to consider 
only a single branch of fermions (for example the right moving one), the single particle 
Green's function is: 

G(X -X',t- t') = ( ^n^^^mkf(x-x') /g^N 

Fourier transforming the equal time Green's function we obtain the average occupations: 

n(fc)^c (f ~^y% n^^g(fc), (32) 
\k — 'mkf\ 

where C is a numerical coefficients and vF^^i^z) is the non-singular part of n{k\ The average 

occupation of the edge excitations must therefore exhibit a power law singularity at mkf 

with exponent m — l. We next present numerical results. From here on we will only consider 

the completely squeezed geometry 7 = 0. 

Fig. 5 shows n(k) for up to 10 electron size systems for u = 1/3. As can be seen, 

there is a remarkable degree of convergence and the shape of n{k) is already clear. We have 
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defined kf{N) for eacli size by liaving exactly A^-states between —kf < k < kj. For each 
size we liave rescaled k in Fig. 5 by kj. Tfie solid curve is n{k) of the Calogero-Sutherland 
jlSl model. It was obtained from a conjecture by Haldane [0 giving the full analytical 
expression of the retarded Green's function for the Calogero-Sutherland model at integer 
coupling (A = m for v = 1/m). Haldane's expression reduces to the equal time Green's 



function recently calculated by Forrester. pO[ Note that in Fig. 5 there appears to be a 



weakly singular pseudo- Fermi surface at kf. The singularity of interest however is not at 
kf but at mkf, i.e. 3kf in this case. Unfortunately, even iV = 10 is not sufficiently large 
to make a direct extrapolation of the exponent possible. That is, we have not approached 
the k = 3kf sufficiently closely (see Fig. 5). However, we extract the exponent as follows: 
we observe that the last possible k for which n{k) is non-zero satisfies {kmax — "^kf) ~ 1/A^ 
(see below). Now if n{k) ^ {k — Skf)"^ then n{kmax) ~ 1/^^, which implies that, in the 
thermodynamic limit the quantity N'^n{kmax) should extrapolate to a finite value. Fig. 6 
shows such an extrapolations clearly demonstrating the {k — Skf)"^ dependence. The value 
of exponent has also been confirmed in the planar disk geometry [pl|]. 

At first sight, the data in Fig. 5 seems to be size independent and to exactly form a 
universal curve. This is not so, there are extremely small finite-size effects. In fact, we have 
empirically obtained the exact expression for the last 4 of n{k) in the tail of the distribution 
as a function of system size: 



14 



^° 2 + 9A^(A^ - 1) 

18(iV- 1) 

n-i - 



n-2 



n^3 



{2 + 9N{N - + 3{N - 2)) 

36{N - l)(9iV2 - SON + 22) 
(2 + 9N{N - 1))(2 + 9(A^ - 1)(A^ - 2))(1 + 3(A^ - 3)) 

4(2 + 3{N - 3))(405iV^ - 2970iV3 + 7A25N^ - 7554Ar + 2698) 



(2 + 9N{N - 1))(2 + 9(iV - 1)(A^ - 2))(2 + 9(A^ - 2){N - 3))(1 + 3(iV - 4)) 

where no = n{kmax), ""--i = n{kmax-i), etc. One might wonder if there is a pattern to be used 
for constructing successively higher terms. While the factors in the denominator appear to 
follow a simple pattern, it is clear that the one polynomial in the numerator rapidly becomes 
complicated and is not easily generalized. However it should be noted that, as ^ oo, the 
ratio of these coefficients become very simple: 



^ = 1 (33) 
no 



^ = 3 (34) 
no 

n-2 



no 

n_3 
no 



6 (35) 
10. (36) 



This sequence makes it clear that these ratios are just the expansion coefficient of (1 — x) ^, 
That is: 

n-p ^ (p + l)(p + 2) 
no 2 

This particular increasing sequence has also been noticed by Wen We note in passing 

that the simple extrapolation procedure in Fig. 6, where we essentially extend the line 
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connecting the points = 10 and = 9, out to = oo, gives the value of 0.22 for the 
coefficient of A^~^, which is remarkably close to the exact 2/9 answer. It should also be 
noted here that one can directly extract the exponent from the analytical expression of hq 
by using the identity 

k — mkf m 

to eliminate in favor of k — 3kf. 

We have also calculated the average occupations for the Boson Laughlin state at z/ = 1/2. 
Here the the exact analytical form of n{k) is known p3 | 



2k f 

n{k) = C'e{2kf - k) In — ^, (37) 

where C is another constant and 6 is the step function. Note that the logarithmic depen- 
dence upon k predicts a linear relation between n{k) and k — 2kf, which is clearly seen in 
Fig. 7. It should also be noted that the increase of n{k) with size at = is highly sugges- 
tive of a weak singularity, consistent with the logarithmic form in Eq. 33. This singularity 
is however an artifact of the squeezed limit and will disappear for non-zero 7 as will the 
singularity at /cj in the Fermi case. 

There are, in the realm of the quantum Hall effect, other states with perhaps more 
interesting edge excitations. We will not address all such states here but confine our attention 
to the pairing ones which have been proposed as candidates for the u = 5/2 Hall state. First 



a spin-singlet state was proposed by us |2^ : 
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"^HR = n(4 - 4f - 4)' n(4 - 4)'Det{-^-^— >, (38) 

i<i i<j i.j {^i Zj ) 

where Det stands for the determinant. 



The other a spin-polarized state was proposed by Moore and Read and later studied 



by Greiter, Wen and Wilczek |£6| 



"^MR = - Zjfm . ^ A , (39) 

i<j \Zi Zj) 

here Pf denotes the Pfafiian. Both the determinant and the Pfaffian act to create pairs 
and are closely related to BCS wave function for singlet and triplet pairing respectively. One 
therefore might expect as in BCS theory a smoothly falling n{k) with no semblance of any 
singularity at kf that was seen for the 1/3 state. This is clearly borne out in Fig. 8 showing 
n{k) for the spin-singlet state for up to 8 electrons. There is indeed a marked difference 
from that of the Laughlin wave function. There is evidently, as before, a reasonable degree 
of convergence already with these small sizes. 

Analogous results for the spin-polarized Pfaffian state [^,|^ is shown in Fig. 9. There is 



a rather sharp drop at kf possibly indicating a weak singularity, i.e. a pseudo- Fermi surface, 

which is probably an artifact of the squeezed limit. In addition, it shows a very slight dip 

at k = 0. Despite their common pairing nature, the two n{k) of Fig. 8 and 9 appear to be 

different. Presumably because of the differences in even (spin-singlet) or odd pairing states 

of the two. A more detailed understanding of these differences will probably have to await 

until a more precise analytical expression for n{k) is obtained. 
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In closing we emphasize again that while for convenience we have studied the edge states 
in the squeezed limit 7 = 0, our results on the edge singularity remain valid for all aspect 
ratios 7 so long as the bulk is an incompressible fluid. As stated earlier only the apparent 
singularities at A; = /c/ for fermions and /c = for bosons are washed out at non-zero 7. 

This work was supported by NSF Grants DMR-9 113876 at CSULA, and DMR-8901985 
at Princeton University. 
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Appendix A 

The first property follows from the antisymmetry of the wave function. We find the 
state with the smallest amplitude by choosing the configuration with lowest power of Zi, 
then lowest power of Z2, etc. We obtain: 

i-f n n ■ ■ ■ ^n^""''^ = i-f n ^t^~'^ (40) 

j=2 j=3 j=2 

From this we immediately pick out the occupation numbers: 

rij — 1 for j — 0, m, 2m, . . . , Nm. (41) 
rij — otherwise. (42) 

Thus, every m^^ orbital is occupied, which is the TT state. Note that TT state in config- 
uration space can be obtained by antisymmetrization of this expression (see below). There 
is another more direct way of seeing this which also will prove point c. We expand the 
polynomial in the following manner: 

- = IK^r -zT + mX,,) (43) 
J{{^T-^T) + rnJ{G,, (44) 

The first term is the TT state with a coefficient of unity. Xij and Gij are appropriate 

remainder terms. It can be seen from the binomial expansion that a factor of m multiplies 

Xij and hence Y[i<j Cij- This proves point c. 

To estabhsh the "squeezability" we write: 
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Af-1 



Eiz^ - = n 



m(N-j) 



i<3 



(45) 



Note the prefactor is related to a similar expression given above (Eq. 7) by the permu- 



tation 



V 



, with the correct sign (— ) , and thus it describes the TT state. 



12 ... N 
N N -1 ... I 

Now consider any pair I with k > I, since in the prefactor zi is raised to a higher power 



than Zk, but in the product Y[i<j{i — Zi/zj), j is greater than i, it then follows that the power 
of zi will increase while the power of Zk will decrease, indicating the pair is being squeezed 



together. 
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FIGURES 

FIG. 1. The charge density plotted versus y the dimension along the length of the cylinder for 
N = 10. Also plotted arc the average occupations of the Landau level orbital placed at the position 
of the guiding centers. The aspect ratios o = are also indicated, in parenthesis, on each plot. 
The series of plots show the evolution of the incompressible fluid from the squeezed cylinder in the 
lower left corner as 7 is increased counter-clockwise. 

FIG. 2. Same as figure 1 except 7 is increased further, stretching the cylinder. The 
charge-density wave state starts to appear at an aspect ratio of 2 and is fully developed at 5. 

FIG. 3. The charge density of the completely squeezed cylinder. The occupation amplitudes 
are now centered at y = since no motion of the Gaussian envelope is possible in the squeezed 
direction. The units of the density in this figure are arbitrary and independent of those in Figs. 1 
and 2. A'^ = 10 in this case as well. 

FIG. 4. An arbitrary energy versus momentum dispersion relation E{k) for a 1-d Fermi gas. 
The Fermi energy is shown by the horizontal line and the Fermi points are labeled. 

FIG. 5. The average occupation of n{k) versus k for the Laughlin state at 1/3 filling for 
up to 10 electron size systems. The /c-axis has been normalized by kf so that k/kf = 2mi/N, 
rrii = —m{N — l)/2, • • • , m{N — l)/2. The solid curve is the exact n{k) for the Calogero- Sutherland 
model. 

FIG. 6. The plot of N'^n{kmax) versus 1/N. An estimate of N'^n{kmax) as A — >■ 00, based on 
the extrapolation of the data, yields 0.22. 

FIG. 7. Same as Fig. 5 but for bosons at u = 1/2. Here there is a hint of the logarithmic 
singularity in the exact result for A^ ^ 00 (see text) at A; = 0. The logarithmic dependence yields 
a linear exponent at 2kf which can be seen in this plot. 
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FIG. 8. Same as Fig. 5 except for the spin-singlet pairing state of Eq. 34. n{k) is similar to the 
single particle occupation in the BCS ground state with no singularity at A;/. 

FIG. 9. Same as Fig. 8, but for the PfafHan pairing state. The diffrence is presumably due to 
S-wave(Fig. 8) versus P-wave pairing. 
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